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£1^, Abstract. The collapsar model for gamma-ray bursts requires three essential ingredients: a massive core, removal 

■ of the hydrogen envelope, and enough angular momentum in the core. We study current massive star evolution 
models of solar metallicity to determine which massive star physics is capable of producing these ingredients. 
In particular, we investigate the role of hydrodynamic and magnetic internal angular momentum transport and 
binary mass and angular momentum transfer. We follow the evolution of rotating single stars and of binary 

| systems that include rotational processes for both stars. Neglecting magnetic fields, we show that the cores of 

massive single stars can maintain a high specific angular momentum (j~10 17 cm 2 s _1 ) when evolved with the 
ly-^ ' assumption that mean molecular weight gradients suppress rotational mixing processes. In binary systems that 

£ — . undergo mass transfer during core hydrogen burning the mass receiving star accretes large amounts of high 

t-H ■ angular momentum material, leading to a spin-up of the core. We find, however, that this merely compensates 

for the tidal angular momentum loss due to spin-orbit coupling, which leads to synchronous rotation before the 
, mass transfer event. Therefore the resulting cores do not rotate faster than in single stars. We show that some 

accreting stars become Wolf-Rayet stars at core helium exhaustion and form CO-cores that are massive enough 
to form a black hole. We also present models that include magnetic fields generated by differential rotation and 
we consider the internal angular momentum transport by magnetic torques. Though magnetic single star models 
O i 1 are known to develop rather slowly rotating cores with specific angular momenta at the end of the evolution 

close to those in observed young pulsars (j~10 14 cm 2 s _1 ), we investigate the capability of magnetic torques to 
efficiently pump angular momentum into the cores of accreting stars. Despite our finding that this mechanism 
works, the magnetic coupling of core and envelope after the accreting star ends core hydrogen burning leads to 
slower rotation (j~10 15 ~ 16 cm 2 s _1 ) than in the non-magnetic case. We conclude that our binary models without 
magnetic fields can reproduce stellar cores with a high enough specific angular momentum (j>3T0 16 cm 2 s _1 ) to 
, ' produce a collapsar and a GRB. If magnetic torques are included, however, GRBs at near solar metallicity need 

, to be produced in rather exotic binary channels, or current dynamo model overestimates the magnetic torques. 

■ But then the problem is that significant angular momentum loss from the iron core either during core collapse or 
' from the proto-neutron star would be required. 
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1. Introduction finally, their afterglows show signatures of the shaping of 

the circumstellar med ium by a massive progenitor star 

From studies of host galaxies of gamma-ray bursts it was J van Marie et al.ll2004lh The most widely used model for 

concluded that they occur in or close to star forming re- GRB production in the context of black hole formation 

gions. Several of them are associated with a very energetic j n a massive sin gle star is the so called collapsar model 

variety of Type Ic supe r nova with broad l ines ("hyper- (|Wooslevlll993all 

nova" , IPrice et all 120021: iHiorth et all l2003j) , thought to 

display the explosion of a massive Wolf-Rayet star. And A collapsar is a massive (M > 35-40 M^. lFrverlll999(l 

rotati ng star whose core collapses to form a b lack 
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(j>3-10 16 cm 2 s- 1 . lMa,cFa,dven fc Wooslevlll 999h an ac- 
cretion disk is formed around the black hole. The accretion 
of the rest of the core at accretion rates up to 0.1 M©s _1 
by the newly-formed black hole is thought to be capable 
of producing a collimated highly relativistic outflow. This 
releases large amounts of energy (~10 51 erg s _1 ) some of 
which is deposited in the low density rotation axis of the 
star. In case the star has no hydrogen envelope, i.e., has a 
light crossing time which is less or comparable to the du- 
ration of the central accretion (about 10s), a GRB accom- 
panied by a Type Ib/c supernova may be produced. The 
collapsar models for gamma-ray bursts thus need three 
essential ingredients: a massive core, loss of the hydrogen 
envelope, and sufficient angular momentum to form an 
accretion disk. 

iHeeer et~ai1 l|2000b|) have calculated models of a 25 M© 
star that could form a black hole by fallback (SN explosion 
occurs while the core forms a neutron star, but so much 
matter fails to escape and falls back onto the neutron star 
that it turns into a black hole). This star ends its life as 
a red supergiant with an iron core of 1.9 Mq, a helium 
core of 8.06 Mq, and a low density envelope of 6.57 M©. 
They found that this star has sufficient angular momen- 
tum to form an accretion disk around the black hole 
which may lead to an asymmetric, jet-driven supernovae, 
or in the case that the star lost its hydrogen envelope, 
even a GRB can result. T he rotating prc-supernova mod- 
els of IHeeer et ail (|2000a|) predicted, for the mass range 
10...20M©, an iron core angular momentum barely suf- 
ficient for the GRB production model through collapsar, 
with a trend of decreasing f inal specific a ngular momen- 
tum for larger initial mass. IHeeer et all (|2000a|) showed 
also that specific angular momentum of the stellar core 
depends significantly on the inhibiting effect of the mean 
molecular weight gradient on rotational mixing processes. 

Only if composition is a efficient barrier for rotational 
mixing and transport of angular momentum, it prevents 
the core from losing most of its angular momentum during 
the evolution. Finally, the rotating pre-supernova models 
of (|Hirschi et all 12004^ predict about twice as much fi- 
nal cor e angular momentum as the models of lHeeer et alJ 
1)20004 . 

These calculations left out the influence of a magnetic 
field which can significantly alter the angular momentum 
transport processes in the star, for example slowing down 
the he lium core of the star during the red supergiant 
phase l|Spruit fc Phinnevll998() . Single star pre-supernova 
models including angular momentum transport b y mag- 
netic t orques, using the improve d dynamo model of | Spruid 
<l2002l). have been p roduced by iHeger et~ ] ll2004air and 
IHeger et all l)2004fj) . with a clear result: although these 
models predict neutron star spins in the range displayed 
by young pulsars, the amount of angular momentum in 
their cores is one to two orders of magnitude less than 
what is required by the collapsar model of GRB produc- 
tion. A possible conclusion the refore might be that, if the 
dynamo model of lSoruitl l(2002l) is qualitatively right, mas- 



sive single stars can not form GRBs within the collapsar 
model. 

Only one in one hundred collapsing massive stars needs 
to produce, however, a GRB in order to explain their fre- 
quency. It seems therefore legitimate to explore collapsars 
in the fra me of massive close bin ary evolution. IWellsteiiil 
(|200l|) and IPetrovic etail l)2004j) showed that during the 
mass transfer phase in a binary system, the secondary (ac- 
creting) star can spin up to close to critical rotation, i.e., 
surface layers of this star can gain large amounts of an- 
gular momentum. This angular momentum can be trans- 
ported inward and increase the rotation rate of the stellar 
core. In this case, one may expect that the effects of mag- 
netic fields would be helpful in spinning up the core, and 
that faster rotating cores might be produced than in the 
non-magnetic models. 

In this paper, we investigate four different types of ro- 
tating massive star models: single and binary models, both 
with and without magnetic fields included. In Sect. we 
briefly explain our numerical methods and physical as- 
sumptions, and in Sect.|3|we explore non-magnetic single 
and binary stars. In Sect.0]wc present our magnetic mod- 
els, first single stars then binary evolution models, and in 
Sect. we summarize our conclusions. 



2. Computational method 

We compute detailed evolutionary rotating models of non- 
magnetic 20 M© and 42 M© single stars, and of a 42 M© 
single star with magnetic field included. Then we calcu- 
late the evolution of a non-magnetic and a magnetic bi- 
nary system with rotating components of Mi.i n =56M© 
and M2 i in=33M© and an orbital period Pi n =6 days. We 
also compute models of a rotating 33 M© star that ac- 
cretes matter at the end of its main sequence evolution, 
for the non-magnetic and magnetic cases. 

W e use the evolutionary code developed by iBraunl 
(|l998l) on the basis of an imp licit hvdrodynamic stellar 
evolution code for single stars jLangerlll99llll998h . 

The treatment of convection a nd se miconvection have 
been described in lLangerl <|l99ll) and [| raun fc Langerl 
l|l995h . Changes in chemical composition are computed us- 
ing a nuclear network including pp chains, the CNO-cycle, 
and the major helium, carbon, neon and oxygen burning: 
reacti ons. More details are give n in IWellstein fc Langerl 
(jl999h and lWellstein et al.l(|200lh . For all models, a mctal- 
licity Z=0.02 is adopted an d we use the OPAL opac- 
ities l lglesias fc Rogers) Il996j) . The abundance ratios of 
the isotopes are chosen to have the solar meteoriti c 
abundance ratios according to iNoels fc Grevessel l)l993|) . 
The change of the orbital period due to mass trans- 
fer and stellar wind mass loss is computed according to 
IPodsiadlowski et alJ l(l992h . with a specific angular mo- 
m entum of the stellar w i nd ma terial calculated according 
to lBrookshaw fc Tavanil l|l993|) . The influence of the cen- 
trifugal fo rce in the rotating models is i mplemented ac- 
cording to lKippenhahn fc Thomas! (|l970() . 



Ste llar wind mass lo s s for stars is calculated accord- 
ing to iKudritzki et alJ l|l989|) . For hydrogen poor stars 
(X s < 0.4) we used a relation based on the empirical mass 
loss ra tes of Wolf-Rayet stars derived by lHamann et alJ 
1)199^ : 

log(M WR / Meyr' 1 ) = -11.95 + 1.5logL/L Q -2MX s .(l) 

Since lHamann fc Koesterkd l)l998j) suggested that these 
mass loss rates may be overestimated, we calculated evo- 
lutionary models in which the mass loss rate given by Eq.^ 
is multiplied by 1/3. 

Also, enhanced mass loss due to rotation is included: 



M/M(v Iot = 0) = 1/(1 - n)< 



(2) 



where £=0.43, n=v Iot /v clit and «^ t =GM(l - T)/R with 
r=i/I/Edd=KL/(47rcGM) is the Eddington factor, G is 
gravitational constant, M is mass, R radius, k opacity, 
Vmt rotating v elocity and i; cr j t critical rotational velocity 
jLangeifllQQsj) . 

The transport of angular momentum in our code is 
formulated as a diffusive process: 
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where v is the turbulent viscosity and i is the specific 
angular momentum of a shell at mass coordinate m. Factor 
(1/2) (dlni/dlnr) vanishes if the gyration constant k = 
i/r 2 does not depend on r, i.e. if all shells are moving 
homologously. 

A parameter / M =0.05 is adopted for sensitivity of the 
rotationally i nduced mixing processes to the M-gradient 
jHeeer et all l2000a|) . As shown by IMevnet fc Maederl 
^129JJh /^-gradients can efficiently suppress rotationally 
induced transport processes. The strenght of this in- 
hibiting effect is descib ed by the parameter / At =0..1 
1 Pinsonneault et "alii 19891) . A value of Z^O.05 reproduces 
the best observations of the enrichement in the surface of 
the main sequence stars with products of CNO process. 

At the surface of the star, the angular momentum con- 
tained in the layers which are lost due to stellar wind gets 
removed from the star: 



J = Mj spc 



(4) 



where j spoc is the average angular momentum at the sur- 
face of the star and M the stellar wind mass loss rate. The 
turbulent viscosity, is determined as the sum of the con- 
vective and semiconvective diffusion coefficients and those 
from rotationally induced instabilities (dynamical shear, 
Solberg-H0iland, secular shear, Goldreich-Schubcrt-Fricke 
instability and Eddington-Sweet circulation). Rotationally 
induced mixing processes and angular momentum trans- 
port through the st ellar interior are described in detail by 
iHeger et all l|2000a|) . 

The code calculates the simultaneous evolution of two 
stellar components of a binary system and computes mass 



transfer within the Roche approximation 

Mass loss from the Roche lobe filling component th rough 

the first Lagrangian point is given bv IRitteil l|l988|) as: 



M = M exp(i?- Ri)/H p 



(5) 



with M =pv a Q I \fe, where e is the base of the natural 
logarithm, H p is the photospheric pressure scale height, 
p is the density, v s the velocity of sound, and Q the 
effective cross-section of th e stream through the first 
Lagra ngian point according to lMever fc Mever-Hofmeisterl 
()l983h . The time scales for synchronization and circular- 
ization of th e binary orbi t as well as spin-orbit coupling 
are given bv lZahnl l|1977|) . The specific angular momen- 
tum of the accreted matter is determined by integrating 
the equ ation of motion of a test particle in the Roche po- 
tential iWenstemll20"0ll accretion stream impacts directly 
on the secondary star). 

We calculate the evolution of binary systems in detail 
until Case AB mass transfer starts. Then we remove the 
hydrogen rich envelope from the primary, until only 5% 
of the hydrogen is left in the envelope. This is the point 
where we assume that the primary shrinks and Case AB 
stops. On the other side, we calculate the secondary star 
assuming accretion with a mass transfer of 10~ 4 M Q yr~ 1 . 
We calculate the Kelvin-Helmholtz time scale of the pri- 
mary on the beginning of Case AB mass transfer: 



t KH = 2 • 10 7 Mi7(£ifln)yr, 



(6) 



where Mi, L\ and R\\ are mass, luminosity and Roche 
radius (in Solar units) of the primary star at the onset of 
Case AB mass transfer. Mass transfer rate is: 



M tr = (Ah - M He )A 



KH 



(7) 



where Mn e is the mass of the helium core with ~5% of 
hydrogen on the surface, i.e. the initial WR mass. We cal- 
culate the orbit change assuming constant mass transfer 
rate calculated by Eq.0and an average accretion efficiency 
(j3) which is the same as for the fast phase of Case A, since 
Case AB also happens on the thermal time scale. Matter 
that has not been accreted on the secondary leaves the sys- 
tem with a specific angular momentum which corresponds 
to th e secondary orbital angular momentum ijKing et alJ 
l200lh . Stellar wind mass loss is neglected. 

We follow the evolution of the binary system, again 
in detail, until the primary finishes carbon core burning. 
We assume that after this the system is disrupted by the 
SN explosion of the primary and we model the secondary 
further as a single st ar. Stellar wind mass loss for red s u- 
pergiants is given bv lNieuwenhuiizen fc de .Tagerl 1 199fl) . 

Magnetic fields ge nerated by dif ferential rotation are 
included according to ISpruitl l)2002j) . The rate at which 
the field is amplified is determined by the differential ro- 
tation. Differentional rotation rate (rotation gradient) is a 
function of a radial coordinate only ('shellular rotation') 
and is given by: 
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The dynamo process requires a minimum rotation gradi- 
ent q m in to operate. 

The effective radial viscosity produced by the magnetic 
field is: 



where 
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r is the radial coordinate, is the rotation rate of the 
star, is the compositional buoyancy frequency, Nt is 
the thermal part of buoyancy frequency, q is the rotational 
gradient and g m i n is minimum rotational gradient neces- 
sary for the dynamo to operate: 



f{l) = 1 - Qmin/q, (q > qmin), 

and 

/(?) = 0, (q < g min )- 



(12) 



(13) 



The factor f(q) causes the stress to vanish smoothly as 
the gradient of the rotation rate approaches the minimum 
value required for dynamo action. 

Fluid motions involved in the dynamo process also im- 
ply a certain amount of mixing. The effective diffusivity 
is given as: 
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where f{q) is defined by Eq.(ll) and 
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3. Models without magnetic field 

3. 1 . Single stars 

It was sh o wn by iMevnet fc Maederl l|l997h and 
Irleerer et al ] i|2000ah that transport of angular mo- 
mentum in stellar interiors depends strongly on the 
efficiency of the inhibition of the rotationally induced 
mixing processes by ^-gradient. When rotationally in- 
duced mixing processes are not inhibited by ^-gradient, 
stellar rotation remains close to rigid rotation. The star 
loses matter from the surface which carries away angular 
momentum. Layers beneath the lost ones expand, and 
due to local angular momentum conservation, spin down. 
The star reestablishes rigid rotation by transporting 
angular momentum from the core to the surface. When 



the rotationally induced mixing processes are inhibited by 
the /i-gradient, angular momentum cannot be efficiently 
transported between the core and the envelope. The 
result is differential rotation between these two stellar 
regions. In this paper we present models which include 
the inhibiting effect of /x-gradient on rotationally induced 
mixing processes. 

We calculate the evolution of a 20 M Q star until the 
end of core carbon burning. This star is a rigidly rotating 
star on the ZAMS (Fig. Q three dots-dashed line) with 
an initial surface velocity of z wf=200 km s -1 , which is 
a typical value for these stars l)Fukudalll982fl . During its 
hydrogen core burning phase, this star loses mass due to 
a stellar winds (~10 _8 Mq yr _1 ). Matter lost from the 
stellar surface carries away angular momentum, the sur- 
face layers spin down, but there is no efficient transport 
of angular momentum between the core and the envelope, 
so the core does not spin down significantly. At the time 
of helium ignition (Fig.^dot-dashed line), the specific an- 
gular momentum of the core is j~10 17 cm 2 s _1 at 3M Q 
(further in the paper, specific angular momentum values 
are always given at 3M Q , since that is the mass of the 
final core before supernova explosion). After core hydro- 
gen exhaustion, the star evolves into a red supcrgiant, its 
core contracts while the envelope expands. This leads to a 
spin-up of the core and a spin-down of the envelope. The 
envelope is convective and rotating almost rigidly with 
a very low rotational velocity (~0.1 km s _1 ). The core 
is rigidly rotating with a maximum rotational velocity of 
~130 km s _1 . The core and the envelope are separated 
by layers that have a large /j,-gradient. This suppresses ro- 
tationally induced mixing and angular momentum is not 
efficiently transported, so that the core is not slowed down 
by the slow rotation of the envelope. At the end of helium 
core burning the specific angular momentum of the core 
is j^5T0 16 cm 2 s" 1 (Fig. fj] dotted line). During further 
evolution, i.e., core carbon burning, the stellar core does 
not lose any significant amounts of angular momentum. 

To investigate the behaviour of higher mass stars, we 
modeled the evolution of a 42 M Q star with the same ini- 
tial surface rotational velocity w sur f=200 km s -1 . The star 
is a rigidly rotating star on the ZAMS with the specific an- 
gular momentum profile shown on Fig.|3(dot-dashed line). 
As we already explained for the previous example, the star 
loses mass due to a stellar wind (^-dO~ 7 M Q yr _1 ). Since 
the stellar wind mass loss rate is one order of magnitude 
higher than for a 20 M Q star, the 42 M Q star loses more 
matter and angular momentum during the core hydro- 
gen burning phase. When helium ignites in the core, the 
specific angular momentum at 3M Q is ^8T0 16 cm 2 s -1 . 
The star becomes a ~38 M Q red supergiant with central 
helium burning in a convective core of ~12 M Q and a con- 
vective envelope of ~20 M©. Due to stellar wind mass loss 
of ~10~ 4 M Q yr -1 , the mass of the envelope rapidly de- 
creases and reaches ~10 M Q at the end of our calculations. 

We follow the evolution of this star until it burned 
~30% of the helium in its core. The specific angular 
momentum of the core (at 3Mq) at that moment is 




Fig. 1. Specific angular momentum profiles of a 20 M Q 
single star on the hydrogen ZAMS (three dots-dashed 
line), when helium ignites in the core, (dot-dashed line), 
when the central helium abundance is 67% (dashed-linc), 
in the moment of core helium exhaustion (dotted line) and 
at core carbon exhaustion (solid line.) 



Fig. 2. Logarithm of the integrated angular momentum, 
J(m)= J Q M j(m')dm' divided by to 5 / 3 , as a function of the 
mass coordinate, to, for a 20 M star for the same evo- 
lutionary stages as shown in Fig. ^ The thin lines give 
a logarithmic scale of levels of constant J labeled with 
log(J/(erg s)). 



~6.5T0 16 cm 2 s _1 . If we assume that the angular mo- 
mentum of the core decreases further during helium core 
burning with the same rate, the specific angular momen- 
tum of the core at the moment of helium exhaustion would 
be ~3.5T0 16 cm 2 s _1 . We see in Fig. ^th a t, there is no 
angular momentum loss from the core during core carbon 
burning, so we can conclude that this star might produce 
a collapsar and in case that the hydrogen envelope is lost 
during the red the supergiant phase, a gamma-ray burst 
can be the result. 

For a rigidly rotating body of constant density po, the 
total angular momentum, ./(to), enclosed below the mass 
coordinate, to, is J( TO )= J q j{m')dm. In Fig.EJand Fig.H 
we plot the logarithm of J(to)/to 5 / 3 as a function of mass 
and a grid of lines of J=const. If angular momentum is 
transported throughout the surface defined by a mass co- 
ordinate to, J(to) drops. If no such transport occurs, J(m) 
remains constant. If a line of constant J is followed from 
one evolutionary stage to another, it can be seen to what 
mass coordinate the angular momentum has been trans- 
ported in the star during the time between the two evolu- 
tionary stages. 



3.2. Binary systems 

A star evolving in a binary system and accreting mat- 
ter from the companion, increases its surface angular mo- 
mentum. If this angular momentum can be transported 
efficiently through the stellar interior, the star may evolve 
into a red supergiant that has a rapidly spinning core with 
sufficient specific angular momentum to produce a collap- 
sar. 

To check if accretion can add enough angular momen- 
tum to the core, we modeled the evolution of a rotat- 
ing binary system with initial masses of -Mi i i n =56M 
and M2 i m=33 M© and an initial orbital period of p; n =6 
days. The binary system quickly synchronizes during the 
main sequence evolution. Due to this synchronization, 
both stars lose angular momentum and their initial sur- 
face rotational velocities are 92 km s^ 1 for the primary 
and 64 km s _1 for the secondary which is much slower 
than the typical va lues for single stars of these masses 
(|Heger et alJl2000a|) . This means that stars in binary sys- 
tems lose a significant amount of angular momentum due 
to synchronization. The angular momentum loss increases 
with the initial orbital period. 




Fig. 3. Specific angular momentum profiles of a 42 M Q 
single star on the hydrogen ZAMS (dot-dashed line), when 
helium ignites in the center (dotted line), and when the 
central helium abundance is G7% (solid line). 



Evolutionary tracks of the primary and the secondary 
star are given in Fig. [3] and Fig. [B] The primary is ini- 
tially the more massive star; it evolves faster and fills its 
Roche lobe during hydrogen core burning. The binary sys- 
tem enters Case A of mass transfer (dotted line, Fig. |5J). 
The primary loses matter with a high mass transfer rate 
(M™ ax ~3.2-10 -3 M Q yr _1 ) and decreases in luminosity. 
During the fast phase of Case A mass transfer, the primary 
loses ~19 M Q and the secondary accretes only about 15% 
of that matter due to the st ellar wind mass loss caused 
by rotation close to break-up (|Petrovic et alJl2004[) . After 
the fast process of Case A mass transfer, the primary 
continues to expand on a nuclear time scale and transfer 
mass onto the secondary star with mass transfer rate of 
^10 -6 Mo yr _1 (slow phase of Case A). At the end of core 
hydrogen burning the primary contracts and thus Roche 
lobe overflow (RLOF) stops. When the primary starts hy- 
drogen shell burning (M=16.6M ) it expands again, fills 
its Roche lobe and Case AB mass transfer starts. We as- 
sume that this mass transfer stops when the WR star has 
5% of hydrogen on its surface (Af WR (A s =0.05)=14.8 M ). 
More details about the evolu tion of this sy s tem u p to 
WR+O phase can be seen in IPetrovic et alJ l|2004h . As 
explained in Sect. [21 we continue following the evolution 



Fig. 4. Logarithm of the integrated angular momentum 
J(m)= J Q M j(m')dm' divided by m 5 / 3 , as a function of the 
mass coordinate m for a 42 M Q star for the same evo- 
lutionary stages as shown on Fig. [3] The thin lines give 
a logarithmic scale of levels of constant J labeled with 
log(J/(erg s)). 



of the binary system with a hydrogen-free WR star. The 
dashed line on Fig. [3] connects the last calculated model 
at the onset of Case AB and the first calculated model 
when the primary is a hydrogen-free WR star. The initial 
hydrogen-free WR star is Mwr=12.8M q with an effec- 
tive temperature of ~1.2-10 5 K, and it loses ~8M Q due 
to WR mass loss during core helium burning. Due to this 
mass loss, the luminosity of the primary decreases (dash- 
dotted line Fig. [5j) . The primary ends carbon core burning 
(dot-dashed line Fig. 0) as ^4.5 solar mass star. The or- 
bital period of the binary system is ~13.5 days when the 
primary explodes in a supernova. 

The secondary begins its evolution as a 33 M core 
hydrogen burning star. It accretes ~3 M during the fast 
phase and ^4 M during the slow phase of Case A mass 
transfer. Due to this mass gain, the luminosity of the 
secondary increases (Fig. EJ. As we already mentioned, 
the primary loses ~1.8M© during Case AB. We assume 
the same average accretion efficiency during this mass 
transfer as during the fast phase of Case A mass trans- 
fer (/3=0.15), since they both take place on the ther- 
mal time scale. This means that the secondary accretes 
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Fig. 5. Evolutionary track of the primary in HR diagram. 
Solid line: core hydrogen burning phase before and after 
Case A mass transfer. Dotted line: Case A mass transfer. 
Dashed line: connection between last calculated model at 
the onset of Case AB and first calculated model when 
the primary is a hydrogen-free WR star. Dash-dotted line: 
core helium burning phase. Dot-dashed line: core carbon 
burning phase. 



~0.25 Mq and after Case AB mass transfer, it is still a 
core hydrogen burning star, more massive than at the be- 
ginning of its main sequence evolution (39.25 Mq). After 
the secondary exhausted all hydrogen in the core, it con- 
tracts and increase its luminosity and effective tempera- 
ture. Due to the increase of temperature in the envelope, 
the star ignites hydrogen in a shell, expands and cools 
down drastically (dash-dotted line Fig. HJ i?«2300 Rq, 
7ofF~3500 K). Meanwhile, the core contracts further and 
the core temperature increases. Helium core burning starts 
when T" c «1.4-10 8 K. We stopped detailed modeling when 
the secondary has ^67% of helium left in the core. 

Fig. shows the evolution of the internal structure 
of the secondary from the ZAMS, through Case A and 
Case AB mass transfer and the red supergiant phase until 
Y c =0.67. The secondary starts its main sequence evolu- 
tion as a 33 Mq star with a convective core of ~20M Q . 
Rotationally induced mixing processes take place in the 
radiative envelope of the star. The secondary then accretes 
~3 M Q during the fast phase and ~4 Mq during the slow 



Fig. 6. Evolutionary track of the secondary in the HR 
diagram. Solid line: core hydrogen burning phase before 
and after Case A and Case AB mass transfer. Dotted line: 
Case A mass transfer. Dash-dotted line: core helium burn- 
ing phase. 



phase of Case A mass transfer. Heavier elements accreted 
on the surface are relocated by thermohaline mixing pro- 
cess and the convective core increases its mass (~25M Q ). 
The secondary becomes a rejuvenated ~39 Mq core hydro- 
gen burning star in a WR+O binary system. The system 
is likely to be disrupted by the explosion of the primary 
while the secondary is still a main sequence star. Details 
of the evolution of the internal structure of the secondary 
after the primary exploded are shown in Fig. |SJ After hy- 
drogen is exhausted in the core, the secondary becomes 
a ~37 Mq red supergiant with central helium burning in 
convective core of ~12 M Q and a convective envelope of 
~20 Mq. Due to stellar wind mass loss of ~10 -4 Mq yr -1 , 
the mass of the envelope rapidly decreases and reaches 
~10Mq at the end of our calculations. 

The modeled binary system starts its evolution with 
both components synchronized to the orbital rotation. 
The surface rotational velocity of the secondary star is 
~64 km s _1 . When the fast Case A mass transfer starts, 
the secondary accretes matter from the primary with high 
mass transfer rates (^10 -3 M Q yr -1 ). This matter carries 
angular momentum and spins up the top layers of the 
secondary star (Fig.^Jl. The mass transfer rate during the 
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Fig. 7. Evolution of the internal structure of a rotat- 
ing 33 M Q secondary from the ZAMS until red super- 
giant phase y c =0.67. Convection is indicated with diag- 
onal hatching, semiconvection with crossed hatching and 
thcrmohalinc mixing with straight crossed hatching. The 
hatched area at the bottom indicates nuclear burning. 
Gray shaded areas represent regions with rotationally in- 
duced mixing (intensity is indicated with different shades, 
the darker the colour, the stronger rotational mixing). The 
topmost solid line corresponds to the surface of the star. 



slow phase of Case A is significantly lower (10 -6 M Q yr -1 ) 
and the surface rotational velocity of the secondary in- 
creases to about 200 km s _1 . For Case AB we assumed 
a mass transfer rate of ~10 _4 Mq yr _1 (Sect. [3} and this 
spins up the surface of the secondary to ^500 km s _1 . 
When the secondary spins up to close to critical rotation 
it loses more mass according to Eq. 2 . High mass loss 
decreases the net accretion efficiency and also removes an- 
gular momentum from the secondary star. The secondary 
star is also spun-down by tidal for ces that tend to syn - 
chronize it with the orbital motion l|Petrovic et aljEoO-^ . 
After Case AB mass transfer, the secondary will synchro- 
nize with the orbital motion in the WR+O binary system. 
After the SN explosion of the primary, the secondary be- 
comes a red supcrgiant with a very slowly rotating enve- 
lope and its surface rotational velocity drops to less than 
0.02 km s _1 . 
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Fig. 8. Evolution of the internal structure of a rotating 
33 Mq secondary from SN explosion of the primary and 
disruption of the system. See Fig. [7| for an explanation of 
the different hatching types. 



Fig. EH shows rotational velocity profiles of the sec- 
ondary star in different evolutionary phases. The surface 
of the secondary gains additional angular momentum dur- 
ing mass transfer, as we already explained. This angular 
momentum is transported through the stellar interior, and 
into the core. The result is that the stellar core spins faster 
on the helium ZAMS (Fig. EH three dots-dashed line) 
than at the beginning of hydrogen core burning (Fig. 1101 
solid line). After this, the core contracts and further in- 
creases its rotational velocity, and the envelope expands 
and slows down. When the central abundance of helium is 
around 67%, the maximum rotational velocity of the core 
is ~100 km s" 1 . 

Fig. 1111 shows specific angular momentum profiles of 
the secondary at different points of its evolution. The spe- 
cific angular momentum of the secondary increases signif- 
icantly due to fast Case A mass transfer (Fig. ^] dot- 
ted line). After this, the secondary loses angular momen- 
tum due to stellar wind mass loss, but also gains certain 
amount through slow Case A and Case AB mass transfer 
(Fig. EI dashed and dot-dashed line). The result is that 
the core has a larger specific angular momentum when 
central helium burning starts than at the beginning of hy- 
drogen core burning. After core hydrogen exhaustion, the 




Fig. 9. Surface rotational velocity of the secondary star. 
Matter transferee! from the primary with high mass trans- 
fer rate ~10 _3 -10 -4 M Q yr _1 during the fast phase of Case 
A and Case AB spins up the surface layers of the sec- 
ondary up to 500-700 km s _1 . During the slow phase of 
Case A mass transfer rate is lower (~1O _6 M0 yr _1 ) and 
surface rotational velocity increases up to —200 km s _1 . 
After Case AB mass transfer the secondary synchronizes 
with the orbital motion in the WR+O binary system. 
After the SN explosion of the primary, the secondary star 
evolves into a red supergiant with slowly rotating envelope 
of -0.02 km s _1 . 

secondary evolves into a red supergiant, the core contracts 
and the envelope expands. This leads to a spin-up of the 
core and a spin-down of the envelope. The specific angu- 
lar momentum of the core at 3M Q is — 5.5T0 16 cm 2 s _1 
(Fig. ^2 three dot-dashed line). The envelope is convective 
and slowly rotating (—0.02 km s _1 ). The core is rigidly ro- 
tating with maximum rotational velocity of —100 km s _1 . 
The core and the envelope are separated by layers that 
have a high //-gradient. Angular momentum is not effi- 
ciently transported through these layers, so the core is 
not slowed down by the slow rotation of the envelope. 
When a third of the central helium supply is exhausted, 
the core (at 3Mq) has a specific angular momentum of 
— 5T0 16 cm 2 s _1 . If we assume that the angular momen- 
tum of the core decreases further during helium core burn- 
ing with the same rate, specific angular momentum of the 
core at the moment of helium exhaustion is expected to 



Fig. 10. Rotational velocity profiles of the secondary star 
on the hydrogen ZAMS (solid line), after the fast (dotted 
line) and the slow (short dashed line) phase of Case A 
mass transfer, after Case AB mass transfer (dash-dotted 
line), when helium ignites in the core (three dots-dashed 
line) and when the central helium abundance is 67% (long 
dashed-linc). 

be — 4-10 16 cm 2 s _1 . As we have already seen for single 
stars, there is no angular momentum loss from the core 
during core carbon burning, so we can conclude that this 
star has enough angular momentum to produce a collapsar 
and, in the case that the hydrogen envelope is lost during 
red supergiant phase, a gamma-ray burst. 

The gain of angular momentum of an accreting star is 
proportional to the amount of accreted matter. We show 
in Fig. El specific angular momentum profiles of the sec- 
ondary (A c =0, T c =1.4-10 8 K) assuming the accretion of 
0.25, 1.4 and 5M Q during Case AB mass transfer. We see 
that the specific angular momentum of the core is higher 
if the star accreted more matter during mass transfer. The 
specific angular momentum of the core at 3 Mq increases 
from 6-10 16 cm 2 s" 1 to 8-10 16 cm 2 s" 1 . 

4. Models with magnetic field 

4.1. Single stars 

We model the evolution of a single 42 M star with an ini- 
tial rotational surface velocity of 200 km s _1 and magnetic 




Fig. 11. Specific angular momentum profiles of the sec- 
ondary star on the hydrogen ZAMS (long dashed line), 
after fast (dotted line) and slow (short dashed line) Case 
A mass transfer, after Case AB mass transfer (dash-dotted 
line), when helium ignites in the core (three dots-dashed 
line) and when the central helium abundance is 67% (solid 
line). 



field included. Fig. 1141 shows specific angular momentum 
profiles in different phases of the evolution. We see that the 
initial (ZAMS) profiles of the star without magnetic field 
(Fig. |21 dot-dashed line) and the star with magnetic field 
(Fig. E| dot-dashed line) arc identical and represent solid 
body rotation. We notice, however, that the magnetic star 
loses more angular momentum from the core during the 
main sequence evolution as well as between hydrogen core 
exhaustion (Fig. El dashed line) and core helium ignition 
(Fig. El solid line). 

The star loses mass by a stellar wind. This removes 
angular momentum from the surface layers and slows 
them down. Since the magnetic torque keeps the star 
close to solid body rotation during the main sequence 
evolution fFig. H5f) . however, angular momentum is trans- 
ported from the stellar interior towards the surface. The 
viscosity due to the magnetic field is a few orders of 
magnitude larger than the one for rotational mixing 
(~10 10-12 cm 2 s" 1 compared with ~10 7-8 cm 2 s _1 ), and 
it can overcome the /i-gradient barrier. Therefore, the sur- 
face of the star is spun up and the core is slowed down. 



Fig. 12. Logarithm of the integrated angular momentum 
J(m)= J Q M j(m')dm' divided by m 5 / 3 , as a function of the 
mass coordinate m for a 33 M Q secondary star for the 
same evolutionary stages as shown in Fig. ^2 The thin 
lines give a logarithmic scale of levels of constant J labeled 
with log(J/(erg s)). 



The surface layers rotate faster than in the corresponding 
non-magnetic star, and the stellar wind mass loss is higher 
(Fig. El- 

Summarizing, the core of a magnetic star loses signifi- 
cantly more angular momentum during the main sequence 
evolution than when magnetic fields are not included. One 
reason is that the /i-gradient does not represent a strong 
barrier for angular momentum transport. And second, be- 
cause the surface layers are spun-up, the stellar wind mass 
loss is larger and this requires more angular momentum 
from the core to keep solid body rotation of the star. 
When hydrogen burning is finished (T C =0.8T0 8 K), the 
specific angular momentum of the core (at 3Mq) is less 
than 3-10 16 cm 2 s" 1 (Fig. El dashed line). 

When the star expands after core hydrogen burning, 
magnetic torques still work toward enforcing solid body 
rotation. The envelope expansion, however, is too fast and 
angular momentum transport is not efficient enough, to 
keep the envelope synchronized with the core (Fig. I15|) . 
Since magnetic viscosity can overcome the ^-gradient, the 
core is slowed down and the envelope is spun up. Stellar 
wind mass loss is enhanced due to faster rotation and is 




Fig. 13. Specific angular momentum profiles of the sec- 
ondary star on the hydrogen ZAMS (solid line) and when 
helium ignites in the core (T c =1.4-10 8 K), assuming the 
accretion of 0.25 M Q (dotted line) 1.4 M© (dashed line) or 
5 M (dot-dashed line) during Case AB mass transfer. 



significantly higher than in the case of the non-magnetic 
star fFig. H6|> . The star loses significant an amount of an- 
gular momentum between hydrogen exhaustion and he- 
lium ignition and the specific angular momentum of the 
core (at 3M©) at the time helium core burning starts is 
only 5-10 15 cm 2 s _1 . 

4.2. Binary systems 

We modeled the evolution of a rotating binary system 
with magnetic fields, initial masses Mi i j n =56 M Q and 
Ai2.in=33 Mq, and an initial orbital period of Pi n —6 days. 
The binary system starts synchronized on the ZAMS with 
initial surface rotational velocities of ^90 km s -1 for the 
primary and ~60 km s -1 for the secondary which is, as we 
already mentioned for non-magnetic models, much slower 
than the typical values for single stars of these masses. 

During the fast phase of Case A mass transfer 
(M tr ~10~ 3 M© yr" 1 ), the primary loses ~20M© and the 
secondary accretes only about 4.5 M© of that matter due 
to the mass loss caused by rotation close to break-up. After 
the fast process of Case A mass transfer, the primary con- 
tinues to expand on a nuclear time scale and to transfer 



Fig. 14. Specific angular momentum profiles of the 42 M Q 
single star with magnetic fields on the hydrogen ZAMS 
(dot-dashed line), when 50% of hydrogen is left in the 
center (dotted line), when hydrogen is exhausted in the 
core (dashed line) and when central helium burning starts 
(solid line). 



mass onto the secondary star, but with much lower mass 
transfer rates (slow Case A, M tr ~10~ 6 M© yr -1 ). At the 
end of core hydrogen burning the primary contracts and 
thus RLOF stops. 

Accretion increases rotational velocity of the surface of 
the secondary star to almost 400 km s _1 during the fast 
phase of Case A mass transfer. The mass transfer rate 
during slow Case A mass transfer is significantly lower 
(10~ 6 Mq yr -1 ), and the surface rotational velocity of the 
secondary increases to about 200 km s _1 . 

Fig. 1201 shows specific angular momentum profiles of 
the secondary at different phases of evolution. The spe- 
cific angular momentum of the secondary increases sig- 
nificantly due to the fast Case A mass transfer (Fig. 1201 
dotted line). Angular momentum is transported more ef- 
ficiently through the stellar interior compared to the non- 
magnetic model, since the incurred magnetic torques are a 
few orders of magnitude more efficient in angular momen- 
tum transport than the rotational instabilities (Fig. I17|) . 
Comparing the specific angular momentum of the non- 
magnetic (Fig. Ill|) and magnetic model (Fig. 120(1 , we no- 
tice that during fast Case A the angular momentum of the 




Fig. 15. Spin period profiles of the 42 M Q single star on 
the ZAMS (solid line), when 50% hydrogen is left in the 
center (dotted line), when hydrogen is exhausted in the 
core (dashed line), and when central helium burning starts 
(dot-dashed line). 



Fig. 16. Stellar wind mass loss of the magnetic (dotted 
line) and non-magnetic (solid line) 42 M star. The star 
with magnetic fields loses more mass because magnetic 
torques transport angular momentum, which spins-up the 
surface layers and enhances mass loss. 



magnetic star increases more than that of the correspond- 
ing non-magnetic star (2-10 17 cm 2 s _1 for magnetic and 
1.25T0 17 cm 2 s _1 for non-magnetic star, ~10 4 yrs after 
fast Case A, at 3M ). 

The accretion stops when the secondary still has al- 
most 50% of the hydrogen to burn in the core. Angular 
momentum is efficiently transported from the stellar core 
to the surface and the /^-gradient can not stop it as in 
the case of the non-magnetic star. During further main 
sequence evolution, the stellar core loses significant an- 
gular momentum and when hydrogen core burning stops, 
the specific angular momentum at 3M Q is 2-10 16 cm 2 s _1 . 
Before helium ignites in the core, the specific angular mo- 
mentum decreases to 6T0 15 cm 2 s . The further evolu- 
tion through the red supergiant phase slows down the core, 
since the envelope slows it down (magnetic viscosity can 
overcome the /i-barrier). 

Apparently, if a star with magnetic fields accretes mat- 
ter half way in the main sequence evolution, it has enough 
time to lose angular momentum before hydrogen is ex- 
hausted in the core. This is why we compute models of 
non-magnetic and magnetic 33 M Q stars that accrete 2 M Q 
at the end of their main sequence evolution (Y" c =0.95). 



We assume an initial surface velocity of these stars of 
~60 km s _1 since they evolve in a binary system. 

When the central helium abundance reaches 95%, the 
core of the 33 M Q star with magnetic fields has a specific 
angular momentum of about 2-10 16 cm 2 s _1 (Fig. [2] dot- 
ted line). On the other hand, the core of a star without 
magnetic fields hardly loses any angular momentum and 
its specific angular momentum is 4-10 16 cm 2 s _1 fFig. 1231 
dotted line). After this, both stars accrete 2 M of matter 
with an accretion rate of 10 _4 M Q yr _1 . We assume that 
in this case stellar wind mass loss is not enhanced by ro- 
tation, so there is no additional angular momentum loss 
due to this effect. Accreted matter adds angular momen- 
tum to the surface layers and this angular momentum is 
then transported through the stellar interior much more 
efficiently in the model with magnetic field. At the end 
of accretion (Fig. 1231 and Fig. [^J dashed line), angular 
momentum has been transported further inward in the 
magnetic star than in the non-magnetic star. If we look at 
the specific angular momentum profiles 10 4 yrs after the 
accretion ended, we can notice that the specific angular 
momentum of the core of the magnetic star increases to 
almost 3T0 16 cm 2 s _1 (at 3M Q ). On the other hand, in 
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Fig. 17. Effective viscosities due to magnetic torques 
(solid line) and due to rotational instabilities (dotted line) 
for a rotating 42 M star at a central hydrogen abundance 
of 50% and at core hydrogen exhaustion. 

the non-magnetic star angular momentum has been trans- 
ported inward but still has not reached the core. 

The magnetic star continues to lose angular momen- 
tum from the core during the remaining main sequence 
evolution (^10 5 yrs), and when hydrogen is exhausted 
in the core the specific angular momentum at 3 M Q is 
1.6T0 16 cm 2 s _1 . As explained above, in the magnetic 
model, the core loses significant amounts of angular mo- 
mentum between core hydrogen exhaustion and core he- 
lium ignition, and when central helium burning starts, the 
specific angular momentum of the core at 3 M Q is only 
6-10 15 cm 2 s^ 1 . 

5. Conclusions 

The single star models presented above confirm that non- 
magnetic stars in the mass range relevant for GRB pro- 
duction through collapsars may retain enough core an- 
gular momentum for a GRB to form. In particular, the 
stars that evolve into red supergiants and become Wolf- 
Rayet stars only at the end of core helium burning can 
avoid a significant core angular momentum loss through 
Wolf-Rayct winds and may retain a specific angular mo- 
mentum of up to 10 17 cm 2 s _1 in their iron cores. In these 
models, an effective core angular momentum loss during 



Fig. 18. Evolution of the internal structure of a rotating 
33 M Q secondary star with magnetic field from the ZAMS 
until red supcrgiant phase F c =0.92. See Fig. [7| for an ex- 
planation of the different hatching types. 



the main sequence is prevented by the suppression of ro- 
tational mixing in regions containing a mean molecular 
weight gradient. Angular momentum loss from the stel- 
lar core during heliu m burning is insigni ficant. Contrary 
to the trend found in iHeeer et all l|2000a|) . we find that a 
42 M Q star may end its life with a larger specific core an- 
gular momentum than a 20 M . The reason is that, for the 
same initial equatorial rotation velocity, the initial specific 
angular momentum in the 42 M Q is significantly larger due 
to its larger ini tial radius. We note that a similar trend has 

been found bv lrlirschi et all l|2004|) . 

The dynamo model of ISpruitl (|2002l) cause a signifi- 
cant angular momentum transport even in the presence 
of mean molecular weight gradients. Ilieger et alJ 1 2004al) 
already found the increased coupling of core and enve- 
lope to result in iron core specific angular momenta of 
the order of 10 14 cm 2 s _1 for stars between 10 and 25 M©. 
Our magnetic 42 M© single star model was only computed 
up to core helium ignition. By then, however, the specific 
core angular momentum wa s already red u ced by a factor 
of 30 from the initial value. iHeger et all l|2004a|) showed 
that during core helium burning the core-envelope cou- 
pling reduces the core angular momentum by another or- 
der of magnitude. The final specific core angular momen- 




Fig. 19. Evolution of the internal structure of a rotating 
33 M Q secondary star with magnetic field from after the 
SN explosion of the primary and disruption of the system. 
See Fig. for an explanation of the different hatching 
types. 



turn in our 42 M star can thus be estimated to fall below 
~10 15 cm 2 s -1 , which will render effects of rotation during 
the core collapse insignificant. 

We also model the evolution of a 56M Q +33M Q bi- 
nary system with an initial orbital period of 6 days. Since 
the binary-enhanced mass loss of the primary leads to ex- 
tremely slow rotation, our attention focuses on the sec- 
ondary star, which is spun-up due to accretion of mass and 
angular momentum. We computed two binary evolution 
sequences, with and without magnetic fields. In both se- 
quences, the initial 33 M Q star grows to about 40 M Q due 
to accretion during its main sequence evolution, and sub- 
sequently evolves into a red supergiant. Before the mass 
transfer, tidal spin-orbit coupling leads to bound rotation 
of both stars and to a specific angular momentum which 
is a factor of 3.. 5 smaller than in a corresponding single 
star. 

The accretion leads to a significant spin-up of the star. 
In the non-magnetic model, the mean molecular weight 
gradient limits the inward diffusion of angular momen- 
tum, and the corresponding spin-up of the core merely 
compensates the tidally induced angular momentum loss. 
The result is a helium core that rotationally decouples 



Fig. 20. Specific angular momentum profiles of the sec- 
ondary star with magnetic fields on the hydrogen ZAMS 
(long dashed line), after fast (dotted line) and slow (short 
dashed line) Case A mass transfer, when all hydrogen is 
exhausted in the center (dot-dashed line), when helium 
ignites (three dots-dashed line) and and when the central 
helium abundance is 92% (solid line). 



from the envelope with roughly the same mass and angu- 
lar momentum as in the case of a 42 M Q single star, i.e., 
with high enough specific angular momentum to produce 
a GRB. 

In the magnetic model, the core spin-up due to accre- 
tion is stronger. It temporarily leads to a core spin rate 
which is factor of 2.. 3 above that of a ZAMS star of com- 
parable mass. Magnetic core-envelope coupling, however, 
has reduced the specific core angular momentum by al- 
most a factor 100 by the time the star has started core 
helium burning. Its final core angular momentum will thus 
be comparable to that of the magnetic 42 M© single star 
model discussed above: probably too small to produce a 
GRB. 

Clearly, the effect of core spin-up due to accretion on 
the final core angular momentum will be larger if the ac- 
cretion occurs later in the evolution of the accretion star, 
as there will be less time to lose again the angular mo- 
mentum gained by accretion. Therefore, we performed the 
numerical experiment of accreting 2 M Q of material on a 
33 Mq at a time when the core hydrogen concentration 




Fig. 21. Evolution of the internal structure of a rotating 
33 M Q secondary star until helium ignition in the core. 
The star accretes 2 M Q when the central helium abun- 
dance is 95%. See Fig.[7|for an explanation of the different 
hatching types. 



Fig. 22. Evolution of the internal structure of a rotating 
33 M Q secondary star with magnetic field until helium ig- 
nition in the core. The star accretes 2 M when the central 
helium abundance is 95%. See Fig. for an explanation of 
the different hatching types. 



was down to 5% (instead of about 40% in the binary evo- 
lution model). The result was again that the core-envelope 
coupling in the magnetic model was overwhelming and 
prevents this scenario to be a realistic option for GRB 
production. 

We conclude that our binary models without magnetic 
field can reproduce stellar cores with a high enough spe- 
cific angular momentum (j>3T0 16 cm 2 s _1 ) to produce a 
collapsar and a GRB. 

If magnetic field is taken into consideration, however, 
GRBs at near solar metallicity need to be produced in 
rather exotic binary channels, or the magnetic effects are 
overestimated in our current models. The first option is 
not implausible, since reverse mass transfer from the orig- 
inal secondar y star onto the prima ry star during its Wolf- 
Rayet phase ( Wellstein et allEoOll) , or late stellar merger 
may lead to an efficient core spin- up. The realization fre- 
quency of such events, however, even though it is difficult 
to estimate, may still be small. The latter would require 
a significant angular momentum loss from the iron core, 
either during collapse or from the proto-neutron star, in 
order to explain the r elatively slow rotati on rates of young 
pulsars (cf., however. iHeger et al.ll2004bl) . 
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